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Abstract 

The control point form (CPF) provides explicit control of physical grid 
shape and grid spacing through the movement of the control points. The 
control point array, called a control net, is a sparse grid-type arrangement 
of locations in physical space with an index for each direction. As an 
algebraic method CPF is efficient and works well with interactive computer 
graphics. A family of menu-driven, interactive grid-generation computer 
codes (TURBO) is being developed by using CPF. Key and unique features of 
Turbol (a TURBO member) are discussed and typical results are presented. 
Turbol runs on any IRIS 4D series workstation. 


1. INTRODUCTION 

The numerical solution of fluid flow problems is dependent on the 
discrete representation of the flow domain. In order to accurately simulate 
complex flows, it is important to use a grid that is boundary fitted, alined 
with the flow, concentrated where flow properties change rapidly, and as 
nearly orthogonal as possible. But geometric constraints often conflict 
with these desired grid qualities. One of the best approaches to a sensible 
compromise between these conflicting demands would be to generate the grid 
interactively by using a method that allows easy and precise local control 
over the grid distribution. 

The control point form (CPF) provides precise control of the physical 
grid shape and grid spacing through the movement of the control points. The 
control point array is a sparse grid-type arrangement of locations in 
physical space with an index for each direction. Each control point has a 
domain of influence on the locally surrounding grid points. This local 
control capability of CPF is invaluable when a compromise is needed in a 
conflict between the desired grid structure and geometric constraints. Two- 
and three-dimensional forms of CPF and earlier demonstrations are presented 
in detail in references 1 to 4, but for completeness and clarity the 
two-dimensional form of CPF is included in this paper. 

A family of menu-driven, interactive grid-generation programs, called 
TURBO, is being developed by using the CPF technique to generate grids for 
numerical flow simulations. This paper focuses on Turbol, a member of the 
TURBO family. Turbol is the very basic grid-generation code of TURBO. 



Other TURBO family members under development are TurboAZ, which will have 
automatic zoning capability; TurboAD, which will adapt the grid to solu- 
tions, grid quality measure, and geometry; and TurboT, which will be spe- 
cialized for turbomachinery. With Turbol a user can interactively change 
the structure of the control net and subsequently change the grid structure. 
Key features of Turbol are discussed with illustrations in subsequent 
sections. Turbol constructs an initial control net (from either a given 
boundary grid or a given field grid that requires modification), has various 
levels of local and global grid control, redefines a block to operate on, 
and performs interactive curve discretization for a boundary grid. Boundary 
curve discretization, if desired, can be performed interactively by choosing 
a curve segment and assigning a desired number of points on the chosen 
segment. A damped bicubic spline (ref. 5) is available for discretization 
of curve segments with sharp corners. TURBO can be used as a stand-alone 
code as well as a supportive and complementary operating tool. 


2. CONTROL POINT FORM OF ALGEBRAIC GRID-GENERATION TECHNIQUE 


The control point array is a sparse grid-type arrangement of locations 
in physical space with an index for each direction. In two dimensions it 
will be denoted by (Cy). As an algebraic method the control point form 
provides explicit control of the physical grid shape and spacing through the 
dynamic movement of the control points. Figure 1 shows an example of a 
two-dimensional control point array. 

A fundamental part of the control point formulation is the construction 
of curves. This construction represents algebraic coordinate generation in 
a single direction between two opposing boundaries (fig. 2). These bound- 
aries are the fixed endpoints of the curves. The remaining control points 
are in the interior of the sequence and are used to control the shape of the 
curve. 

From any one control point the two neighboring points on each side form 
local segments within the entire piecewise linear curve. The local segments 
about each control point define a change in direction as shown in figures 1 
and 2. The rate of change is determined by the two linear connections 
attached to the given control point as the curve assumes the respective 
directions between each pair of control points. 

A continuous direction field is obtained in a smooth manner by inter- 
polation. The independent variable for the interpolation is the curve para- 
meterization. The interpolated result defines the field of vectors that are 
tangent to the desired curve. Simply stated it is an interpolation of first 
parametric derivatives. This determines a smooth first derivative of the 



2 




entire curve. The desired curve is then obtained by a parametric integra- 
tion. The integration here is taken so that the curve connects the speci- 
fied endpoints. 

To state the result mathematically, let C lj ' C 2j ' * ’ ’ ' C Mj be the given 
sequence of M points in space as shown in figure 2 (e.g., C 12 ,C 22 , . • • ,C 52 for 
j=2 in fig. 1); let r be the curve parameterization; let Ej(r) be the 
position at r along the desired curve; let r k , . . .r 2 , . . . ,r M1 be the 
successive parametric locations used to interpolate the directions of (C 2j - - 

C lj)'( C 3j " C 2j )'•■•< ( C Mj " C (M-l)j)' and 1et l^i < 1^2 ' ’ ’ ‘ ' ifyi-l be the cor - 
responding interpolation functions that separate each successive pair of 

segments by assuming a nonzero value at the associated location while 
vanishing at the remaining locations for interpolation. In two dimensions 
C . = (x .,y .) and E,(r) = [x,(r), y H (r)]. With the notation the desired 

uj UJ UJ \J J J 

curve is given by 

N-l 

(r) = C M + £ (r > * C'j] for j = 1,2 M U) 


where 


s„(r) -f^OOd* (2) 

Local interpolation functions are used to obtain local grid controls. 
With local functions the movement of a control point results in an alter- 
ation of the constructed curve that is restricted to a local region about 
the point. The remaining regions are unaltered. The simplest local inter- 
polants are the piecewise linear functions that do not vanish over at most 
two intervals defined by Tj < r, <...< r M1 . With the uniform partition r k = 
k for simplicity, an explicit form of the normalized interpol ants , which 
yields G k (r N _ 1 ) = 1 for all k, can be given by 


?i(r) 


2(2 - r) for 1 < r < 2 
0 for 2 i r ^ (N-l) 


iM r ) 


0 for 1 < r < (k - 1) 

(r - k) +1 for (k - 1) < r < k 
(k - r) +1 for k < r < (k + 1) 

0 for (k + 1) £ r < (N - 1) 


0 for 1 < r < (N - 2) 

2(r - N + 2) for (N - 2) <[ r < (N-l) 


(3) 


(4) 


( 5 ) 
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Similarly, a 
index value of i : 

F^t) 


curve can 





be constructed for the other direction for the 
H p (t) (C i>ptl - C ip ) for i = 1,2, ... ,N ( f 


where t 
interpolants. 


is the curve parameterization and H«(t) is the integration of the 

The tensor product form depends only upon C^j and is given by 


T(r, t) = e x (r) 


+ £ 


F ptl {r) - Kfir)] 


(7) 


or alternatively by 


T(r, t) - ^ < t) + £ -*;(<:>] ( 8 ) 

«-l 

These two expressions are equivalent. The tensor product matches Ej or F i at 
the extremities of i and j (e.g., at corner points in fig. 1). 

When boundaries are to be specified, the corresponding data appear at 
the extremities of the values for r and t. Since the coordinate trans- 
formations are generally expressed in the form of a vector for the desired 
positions of all points in physical space, it is convenient to express the 
boundary specifications in terms of the position vector. Thus, the bound- 
aries are denoted by P(l,t), P ( N- 1 , t ) , P(r,l), and P(r,M-l). In order to 
include the boundaries, the multisurface transformation is performed again 
as before, but now with the actual boundaries inserted. This results in a 
modification of T for the r and t directions, respectively. In each such 

directional construction the actual boundaries appear as end conditions for 
the corresponding variable and the remaining boundaries are solely generated 
by the control points. Thus, by subtracting T from the sum of both direc- 
tional constructions, the actual boundaries become end conditions for each 
variable. This process follows a Boolean sum format and upon simplification 
becomes 


Q(r , t) = T(r,t) + a t [l - Gj (r)] [P(l , t) - F L (t)] 

+ a 2 G N _!(r)[P(N - l,t) - F N (t)J 
+ a 3 [l - H i ( t ) ] [P(r , 1) - Mr)] 
+ Q 4 Hm_i (t)[P(r ,M - 1) - E„(r)] 


where each of the four terms following the tensor product T(r,t) represents 
a transfinite conformity to a boundary when each boundary on-off switch is 
1. When any a- is set to 0, the corresponding boundary becomes available 
for free-form modeling by means of the control points. In the order of 
appearance the boundaries are for r = 1, r = N - 1, t - 1, and t = M - 1. 
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3. IMPLEMENTATION OF CPF (PROGRAM Turbol) 


3.1 Construction of Initial Control Net 

Turbol first constructs an initial control net from either given geom- 
etry data, a given boundary grid, or any given field grid that one wants to 
modify. Once a control net is constructed, the grid control features of 
Turbol are ready to operate. When geometry data are given, they are usually 
not sufficient to be used as a boundary grid, as for example in figure 3(a). 
In such a case Turbol interactively chooses a boundary curve segment, 
assigns the desired number of grid points, and stretches them as needed, as 
shown in figure 3(b). An initial control net is then constructed by an 
interpolation (fig. 4). Then from the boundary grid and the initial control 
net, the initial grid is calculated as shown in figure 5. If the field grid 
in figure 5 is given instead of the geometry data or the boundary grid, an 
initial contol net similar to the one in figure 4 is constructed by "attach- 
ment* (ref. 3) to the given grid. Then Turbol produces a field grid, which 
is similar to the given one, from the boundary grid of the given data and 
the control net. 



(a) Given geometry data. (b) Boundary grid. 


Figure 3.-Boundary grid generation. 





(a) Entire control net. (b) Zoomed front and rear portions. 

Figure 4.-Construction of initial control net by interpolation. 
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(a) Entire grid. (b) Zoomed front and rear portions 


Figure 5.-Calculation of initial grid. 

3.2 Primary Turbol Features 

Using workstation mouse keys, users can select any control point in the 
control net. It can then be translated to a new position by translating the 
mouse. Each control point has a domain of influence on the locally sur- 
rounding grid points. The change of a control point location does not 
affect the grid points outside the domain of its influence. The point-by- 
point control of the control net is the basic control feature of Turbol. 

This basic control capability of CPF in two and three dimensions is illus- 
trated in references 2 and 3, respectively. 

For more precise control one needs to use more control points. Control 
points can be added or subtracted from a control net interactively by 
selecting the menu item "Change Number of CPs" and then using the work- 
station mouse as instructed on the screen. The number of control points 
between the bottom and top boundaries was increased from five in the initial 
control net to seven in the control net of figure 6. 

The point-by-point control of a control net can be a very tedious 
process. In order improve the efficiency of the interactive process and to 
give convenience to users, global control capabilities have been added to 
Turbol. One example of this is making the control net orthogonal to the 
boundaries. One only needs to choose menu items. Then Turbol finds normal 
vectors to the boundary and makes the control net orthogonal to the bound- 
ary. The control net shown in figure 6 was changed to the new one shown in 
figure 7 by choosing the menu item "Make Control Net Orthgnl to Bndry." The 
new control net is orthogal to the boundary. This modification in the 
control net structure changed the grid structure from the initial one in 
figure 5 to the new grid, which is orthogonal to boundaries as shown in 
figure 8(a) . 




Figure 6.-Controi net with increased control Figure 7.-Control net orthogonal to boundary, 
points. 
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(a) Orthogonal grid before stretching. (b) Orthogonal grid after stretching. 

Figure 8. -Grid stretching. 

Another global control capability is for grid stretching. Grid 
stretching can be accomplished with Turbol in two ways: through the 

stretching of a control net or through direct grid stretching. Direct grid 
stretching using a stretching function results in more satisfying results, 
especially for viscous flow simulations. The degree of stretching can be 
changed interactively with Turbol. When the grid is stretched as shown in 
figure 8 (b), Turbol automatically constructs a new control net by reattach- 
ment. Further refinement of the grid structure is possible. As an exam- 
ple, the number of control points was increased further from seven to nine 
between the bottom and top boundaries. The new control net was then made 
orthogonal to the boundary. This time the orthogonality was extended to the 
second control points from the boundary, as shown in figure 9. A new grid 
was calculated and is shown in figure 10. An entire view, a zoomed front 
portion, and a zoomed rear portion of the new modified grid are shown in 
parts (a), (b) , and (c) of the figure, respectively. 

Turbol can choose a control point, C^, as a hinge point and choose 
another control point, C^, as a current control point (an intersection of 

two highlighted control net lines on the workstation screen) and can move 
several control points (n - 1 points) between the two simultaneously as if 
they are on a stretched rubberband between and Cy +n . This feature is 

particularly useful when this operation is needed for several values of i as 
illustrated in reference 2 . 

3.3 Geometry Manipulation Using Free-Form Boundary 

The default boundary condition of Turbol is a Dirichlet boundary con- 
dition (i.e., the initial boundary grid distribution is fixed). But if 
desired, the boundary grid can be allowed to move to meet the Neumann bound- 
ary condition that is specified by the orthogonal control net, for instance, 
or to change the geometry itself. To allow the boundary grid points to 
move, select the menu item "Free Form Boundary," which is equivalent to 
turning off one of the boundary on-off switches (a^) of equation (9). For 
illustration, eight control points (Cj, 4 and C 12 i...C 12 4 ) °f the 

control net in figure 9 were moved to the new positions £nown iri fig- 
ure 11(a) by using the point-by-point control feature. The menu item "Free 
Form Boundary was chosen to set a, = 0 and a new grid was calculated. 

Figure 11(b) shows the change of the bottom boundary due to the use of the 
free-form boundary and the domain of influence of the eight control points 
whose positions were changed. 
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(a) Entire control net. 



(b) Zoomed front and rear portions. 
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(a) Entire grid. 




(b) Zoomed front portion. (c) Zoomed rear portion. 

Figure 10 -Grid for extended orthogonality. 




(a) New control point positions. (b) Change in bottom boundary. 

Figure 1 1 -Manipulation using "Free Form Boundary." 
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(a) Before subblock operation. 
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(b) Modified subblock control net. 



(c) After subblock operation. 


Figure 12,-Subblock operation. 


3.6 Input and Output 

Input data for Turbol can be geometry data, a boundary grid, or a 
field grid for attachment. Turbol can save grid and control net data in a 
file at any stage of an interactive session. These saved data can be used 
for flow simulation or can be read in for a future interactive session. 


4. CONCLUDING REMARKS 

Turbol demonstrates its interactive grid-generation and grid mani- 
pulation capability when the control point form is used. Turbol can be used 
as a stand-alone code as well as a complementary operating tool for explicit 
grid manipulation. 
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